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SUMMARY 

An experimental Investigation of the three-dimensional flow field through 
a low aspect ratio, transonic, axial-flow fan rotor has been conducted using an 
advanced laser anemometer (LA) system. Laser veloclmeter measurements of the 
rotor flow field at the design operating speed and over a range of through-flow 
conditions are compared to analytical solutions. The numerical technique used 
herein yields the solution to the full, three-dimensional, unsteady Euler equa- 
tions using an explicit time-marching, finite volume approach. The numerical 
analysis, when coupled with a simplified boundary layer calculation, generally 
yields good agreement with the experimental data. The test rotor has an aspect 
ratio of 1.56, a design total pressure ratio of 1.629 and a tip relative Mach 
number of 1.38. The high spatial resolution of the LA data matrix (9 radial by 
30 axial by 50 blade-to-blade) permits details of the transonic flow field such 
as shock location, turning distribution and blade loading levels to be Investi- 
gated and compared to analytical results. 


INTRODUCTION 

Advancement In the aerodynamic design of gas turbine engine fans, com- 
pressors and turbines have historically resulted primarily from an empirical 
approach. Although this approach has led to more efficient, lighter-weight 
engines, major new advancements In component performance will most likely not 
result from empiricism alone. A fundamental understanding of the complex flow 
phenomena which occur within advanced turbomachines Is essential If significant 
advancements In performance are desired. As the cost of building and testing 
gas turbine engine components Increases so does the need to reduce the number 
of redesigns necessary to achieve the desired performance. The need for 
redesigns Is Inherent In design systems which rely heavily on empirical corre- 
lations. There exists a need to develop design systems which rely less on 
empirical correlations and more on procedures which are capable of calculating 
those aspects of the Internal flow field for which the correlations were devel- 
oped. The development of these new design systems will require extensive ex- 
perimental Investigations of the Intra-blade flow field to determine their 
validity. 



In addressing this need, the NASA Lewis Research Center, In cooperation 
with the Army Propulsion Laboratory, has undertaken a program designed to In- 
crease the basic understanding of the complex Internal flows of gas turbine 
fans and compressors. The program relies primarily on the use of nonlntruslve 
optical techniques to evaluate the Intra-blade flow field of transonic compres- 
sor blading. Concurrently, advanced computational techniques are being devel- 
oped which Incorporate the three-dimensional flow field effects of this class 
of turbomachine. A combination of experimental observation and computational 
experience should lead to a better understanding of these complex flows and 
ultimately to major strides In gas turbine design technology. 

Previous researchers have presented comparisons between three-dimensional 
flow calculations and Internal flow field measurements (refs. 1 and 2). They 
have generally found that Invlscld codes model the Inlet to transonic rotors 
quite well but tend to over predict the deceleration through the rotor. This 
Is usually attributed to lack of the streamline displacement effects due to 
boundary layer growth. Those who have added boundary layer calculations to 
three-dimensional analyses (refs. 3 and 4) have found that the flow field Is 
substantially altered by addition of the displacement thickness to the solid 
surfaces. The technique used In reference 3 Included simple models for shock- 
boundary layer Interaction, laminar-turbulent transition, and separation 
effects. The boundary layer calculations were performed only on the blade 
surfaces and yielded better agreement with experimental data than did the cal- 
culations without boundary layer effects. In reference 4 results were shown 
for the same rotor analyzed In reference 1 and demonstrated the extent to which 
Inclusion of a simple boundary layer model on the blades and endwalls could 
alter the flow field In a transonic rotor. 

The focus of the present work Is the Investigation of a detailed set of 
optical measurements of an advanced transonic, axial-flow compressor flow field 
obtained under the current NASA/Army research program. These data, coupled 
with the results of an Invlscld, three-dimensional flow field calculation, 
provide the basis on which the Internal flow features and trends are Inter- 
preted herein. Strazlsar (ref. 5) provides a more fundamental overview of the 
rotor flow physics, Including the three-dimensional nature of the shock fronts 
based on these same data. 


COMPRESSOR ROTOR 

The research vehicle used In the present work was the first-stage rotor 
of a NASA Lewis designed two-stage fan (ref. 6). The first-stage rotor was 
designed as a 22 blade, low aspect ratio (1.56), damperless replacement for 
the original 43 blade, high aspect ratio (2.94), dampered rotor (ref. 7). The 
new blading met or exceeded Its design total pressure ratio and adiabatic effi- 
ciency goals of 1.629 and 0.896, respectively, when tested In the two-stage 
configuration. The Inlet relative Mach number at the rotor tip was 1.38 at the 
design tip speed of 428.9 m/s and design mass flow rate of 33.25 kg/s. 

The rotor was tested In the present work without Inlet guide vanes or 
downstream stators so that blade row Interactions would not be present. 

Figure 1 shows the measured rotor performance characteristic for a constant 
design speed condition. Also shown In the figure, for reference purposes; are 
data from the redesigned two-stage fan Investigation (ref. 6). Herein, all 
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future reference to rotor performance will deal solely with the Isolated rotor 
data. 


The entire experimental research program conducted with the Isolated rotor 
configuration of this low aspect ratio fan covered a range of data points from 
max flow (MF) to near stall (NS); however, the most detailed Information was 
obtained at the peak efficiency (PE) and near stall points. Most of the com- 
parisons between experiment and computation will be restricted to these two 
flow points. 


INSTRUMENTATION 

A fringe laser anemometer system (LA) was the principal Instrument used 
during the experimental evaluation program. This Instrument Is a single- 
channel, dual beam system with an on-axis backscatter collection scheme and 
has been previously described In detail (refs. 8 and 9). Data aqulsltlon and 
storage was accomplished via a dedicated minicomputer while post processing 
and graphical output was handled by a large central computer. Access to the 
compressor flow path was through a 3 mm thick glass window which extended far 
enough foreward and aft of the test rotor to enable free-stream and blade wake 
measurements. The window curvature closely conformed to the rotor outer flow 
path, thereby minimizing disturbances to the tip region flow. Fluorescent 
seed particles, whose nominal diameter was 1 to 1.4 ym, were spray atomized 
Into the flow stream through a 6 mm diameter tube located over 60 cm upstream 
of the rotor leading edge. In addition to LA measurements, conventional In- 
strumentation was used to record compressor overall pressure and temperature 
levels. This Information was used to calculate the compressor overall perform- 
ance characteristic as well as to control the on-line operating condition set 
point. 


DATA COLLECTION 

The precise rotor-relative location of each LA velocity measurement was 
determined using a pulse from a variable frequency oscillator which was phase- 
locked to the rotor rotational speed. All of the data presented herein were 
collected with the oscillator set at 50 pulses per blade passage and a maximum 
blade survey range of 17 passages. The blade geometry was used to set the 
first measurement position (shaft position) to coincide with the suction sur- 
face of a known blade. Since the data collection windows were equally spaced 
with no circumferential gaps, some data were acquired at shaft positions which 
were physically within the rotor blade. Data collected In these areas were 
handled In the off-line reduction of the data, to be discussed later. This 
method of surveying across a number of blade passages allowed for a larger 
overall data collection window as well as provided the ability to Investigate 
passage to passage flow variations. It Is not within the scope of the present 
work to look at such passage dependency among the data, and thus all of the 
blade passages were ensemble averaged to yield an "average" blade passage data 
set. A typical data run consisted of collecting 60 000 velocity measurements 
within the entire data collection window at each axial survey location. This 
yielded approximately 70 measurements at each of the 850 shaft positions or 
1200 averages for the 50 "average" shaft positions. Time to collect the 
60 000 measurements was approximately 60 sec for each of the two velocity com- 
ponents measured at each shaft position. The two angular orientations of the 
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LA probe volume used to obtain the velocity components was predetermined from 
an Initial screening run at each LA survey span location and flow rate. The 
pltchwlse-averaged ,f low angle at several axial stations along the survey 
streamsurface was obtained from this screening run. The two LA fringe orien- 
tations used In the data runs were derived by bracketing this mean flow angle 
by twenty degrees. 

The laser anemometer survey locations along with the meridional view of 
the computational grid used by the three-dimensional analysis code are shown 
In figure 2(a). The LA blade-to-blade matrix Is shown In figure 2(b). LA 
survey locations, denoted by solid circles In figure 2(a), were axially dis- 
tributed at specific locations based on the axial projection of the local chord 
length along design streamlines to provide data coverage at 20 percent Inter- 
vals upstream, 2.5 percent Intervals around the leading edge, 10 percent Inter- 
vals within the blade and five evenly spaced locations between the rotor 
trailing edge and the position of the stage configuration stator leading edge. 
Survey points were also Included at the conventional Instrumentation survey 
planes (stations. 1 and 2), as shown In the figure. The radial distribution of 
survey points was along design streamlines which pass geometrically through 
desired percent span positions at the rotor trailing edge plane. This radial 
position designates the particular LA survey streamline such as 10, 20, 30 
percent and so on. The nine radial Immersions used In the experimental Inves- 
tigation were spaced at 10 percent span Intervals, beginning at 10 percent 
from the rotor tip and coincided with nine of the three-dimensional computa- 
tional planes. 


NUMERICAL ANALYSIS 

The numerical analysis of the rotor was conducted using a three- 
dimensional computational procedure developed by Denton (ref. 10). This pro- 
cedure was based upon a finite volume technique. The code was modified 
(ref. 4) to Include the effects of boundary layer displacement using an Injec- 
ted mass technique, which allowed flow to pass through the walls In proportion 
to the rate of growth of the boundary layer displacement thickness. The bound- 
ary layers were calculated as If the grid lines In the streamwlse direction 
were actual flow field streamlines. Boundary layers were calculated on both 
the blade surfaces and on the endwalls which were assumed to rotate with the 
blading. The displacement thickness was computed from the von Karman Integral 
boundary layer equation using a constant shape factor (1.5) and constant skin 
friction coefficient (0.005). Setting the shape factor to a constant value 
limits the growth of the displacement thickness In regions of sharp adverse 
pressure gradients and could significantly effect the results. However, In- 
creasing the complexity of the computational procedure could lead to signifi- 
cant Increases In computer time. Thus, the decision was made to evaluate the 
present rotor using the simplified procedure In order to assess Its effect on 
the flow field solution. 

The upstream boundary was located 1.5 times the tip axial chord upstream 
for the tip and 0.7 times the hub axial chord upstream for the hub. At the 
upstream boundary the total pressure, total temperature, absolute tangential 
velocity, and meridional flow angle were held constant at each spanwlse grid 
line and assumed to be uniform In the tangential direction. The only down- 
stream condition specified was the static pressure at the hub, and It was 
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assumed to be constant In the tangential direction. The static pressure grad- 
ient from hub to tip at the downstream boundary was calculated assuming no 
curvature of the streamlines In the hub-to-shroud plane which was appropriate 
for this case based upon the flow path geometry. The downstream boundary was 
located twice the tip axial chord downstream at the tip and one hub axial chord 
downstream at the hub. 

All of the solutions presented were obtained with a computational grid of 
21 points blade-to-blade, 95 points In the streamwlse direction with 50 points 
from blade leading to trailing edge (fig. 3) and 11 points from hub-to-shroud. 

A large number of computational points was used In the present work so that 
any discrepancies between analysis and data could not be attributed to use of a 
coarse grid. The peak efficiency point was repeated with the number of blade- 
to-blade points reduced from 21 to 11. The solution showed that the major dif- 
ference between the coarse grid solution and the fine grid solution was a small 
change In location of the shock and In the relative Mach number ahead of the 
shock at the transonic sections. The shock also tended to be smeared slightly 
more with the coarse grid than with the fine grid. Evaluation of the computa- 
tional method In terms of sensitivity of the code to grid size Is not within 
the scope of the present work and Is mentioned only to provide the reader with 
an Indication of Its Importance for this rotor at this one flow point. The 
sensitivity of the solution to blade-to-blade grid spacing was not checked for 
the other flow conditions analyzed. Streamwlse grid spacing was constant on 
the blading but was expanded upstream and downstream In order to Isolate the 
boundaries from the leading and trailing edges with as few points as necessary 
In an attempt to capture the correct upstream wave strength and to limit the 
effect of a uniform circumferential pressure at the downstream boundary. In 
all cases multigrid (7 by 12 by 3) with a nonuniform tlmestep was used. 

The boundary layer displacement thickness was updated every 50 tlmesteps 
with 50 percent of the new displacement thickness added to 50 percent of the 
previous value. Fractional updating of the boundary layer was necessary since 
full correction tended to produce Instability. 

Convergence was assumed to occur when all the following conditions were 

met: 


(1) The maximum change In axial velocity anywhere In the flow field 
divided by the root-mean-square of all the absolute velocities In the flow 
field did not exceed 0.01 percent. 

(2) Mass flow at each streamwlse grid line agreed with the Inlet mass 
flow to within 0.3 percent. 

(3) Inlet mass flow did not change by more than 0.3 percent over several 
hundred tlmesteps. 

(4) Mach numbers at four selected positions on the rotor suction surface 
did not change by more than 0.02 over several hundred tlmesteps. 

The foregoing conditions listed above are rather stringent for this case 
and may not be required for other blade rows. The 0.01 percent restriction on 
axial velocity change could be relaxed to 0.02 In this case with little change 
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In the solutions. The criterion for change In Inlet mass flow rate was neces- 
sary for this blade row operating at near Its peak efficiency point. For oper- 
ation at lower than peak efficiency flow, small changes In mass flow rate are 
not as Important and that condition could have been relaxed. During Initial 
runs to establish the proper pressure ratio In order to set the peak efficiency 
flow, several solutions were obtained at flows very near that for peak effi- 
ciency (mass flow ratios of 0.984 and 0.998 versus the peak efficiency value 
of 0.992). At 10 percent span from the tip for both flows the shock on the 
suction surface moved approximately 15 percent of axial chord from the position 
obtained at peak efficiency. For the higher flow, the shock on the pressure 
surface moved rearward by about 20 percent of axial chord and Increased In 
face Mach number from about 1.23 to 1.36. 

The code as described In reference 4 did not satisfy the boundary condi- 
tion at a solid surface requiring the velocity normal to the surface to be 
zero. Because of this, there are areas In the flow where the velocity does 
not have a zero normal component to a solid surface. This Is particularly 
true around the leading edge of a blade and at points where a strong shock 
hits the solid surface. All solutions presented herein were produced with 
this condition at the solid surfaces. In order to check the effects on the 
solutions presented, the normal velocities were set at zero (ref. *11) at the 
solid surfaces and the flow point at peak efficiency repeated. Very little 
difference was found In the final solutions when judged In terms of line plots 
of near-blade relative Mach number and flow angles. The flow angle did conform 
to the blade angle at the surface but at one grid point off the blade the flow 
angles for the two cases with and without the zero normal velocity condition 
Imposed were very similar. The major differences occurred at the percent span 
locations near the hub, where the leading edge Is thick. Mach numbers on the 
suction surface were slightly higher than those obtained without the normal 
velocity set to zero. The shock location at all spanwlse stations on the suc- 
tion surface was the same for both cases. The numerical error In total pres- 
sure Increased near the hub sections when the zero normal velocity condition 
was Imposed. 

For the computations with the boundary layers Included, tlmesteps to con- 
vergence ranged from 2486 for the peak efficiency case (where the location of 
the shock Is very sensitive to even small mass flow changes) to 1300 'at the 
maximum flow condition. Computation times on a Cray-1 computer were 17 and 
9 min, respectively. Computation time on an IBM 370-3033 would be approxi- 
mately 10 times as long. No appreciable Increase In time resulted from use of 
the boundary layer calculations. The solution actually converged somewhat 
better with the boundary layer calculations than without. 


DATA PRESENTATION 

Due to the large quantity of experimental data collected, over 32.5 
million velocity measurements per operating point, data presentation represents 
a major effort and concern. To compare the overall flow features observed In 
the experimental and analytical data, contour maps of relative Mach number are 
most useful. Experimental data, however, do not usually produce good quality 
contour maps. The background noise, nonuniform survey location spacing, prob- 
lems associated with LA blade flash and data drop-out all contribute to poor 
quality plots. Analytical data, on the other hand, has very few of these prob- 
lems as evidenced by the quality of such contour maps. 
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To avoid constructing contours of experimental data by hand for each flow 
condition surveyed, a method has been developed to "enhance" the LA data, with- 
out altering Its Informational content, such that It becomes suitable for com- 
puter generated contour plots. The following procedure was used to enhance the 
data before presentation: 

(1) Each blade-to-blade data window was first screened to confirm that 
there were at least one hundred LA measurements obtained and that the window 
circumferential extent did not overlap the pressure surface of a blade. If 
either condition was not satisfied the recorded velocity value was set equal 
to zero. 

(2) An Interpolation procedure was then used to establish values to fill 
the zero velocity data windows using Information from adjacent survey points 
In both the streamwlse and pltchwlse directions. This method allowed Informa- 
tion from both directions to Influence the final Interpolated value. Once the 
areas of zero velocity had been Interpolated, the pltchwlse data was smoothed 
such that a maximum change due to the smoothing was limited to 0.5 percent of 
the original value. This amount of smoothing removed the high frequency noise 
In the data without distorting the highly desirable velocity gradient regions 
of the flow field. 

(3) These nonzero velocity data were then used to Increase the stream- 
wise data density by Interpolating along preferred directions. Streamwlse 
data density of 2.5 percent of axial chord was necessary to achieve the desired 
contour map quality. 

,0 

(4) Finally, these data were smoothed In the chordwlse direction to 
eliminate small oscillations due to acquisition noise In the original data as 
well as possible peaks due to the Interpolation procedure. The smoothing pro- 
cedure permits a maximum change at any point of 0.5 percent of the unsmoothed 
value. 

Figure 4 shows contours of relative Mach number for the 30 percent span 
Immersion survey location at peak efficiency operation. Figure 4(a) represents 
a relative Mach number contour map using only the raw LA data Including zero 
velocity values and blade Interior values. Figure 4(b) shows the Improvement 
that can be achieved with simple removal and replacement of "bad" data and a 
small amount of smoothing. Figure 4(c) shows the effect of Increasing chord- 
wise data density where data has been Interpolated to every 2.5 percent axial 
chord. Finally, figure 4(d) shows the contour map which results from the 
screened, chordwlse Interpolated, smoothed data. As shown, the enhancement 
procedure enables reasonable contour maps to be constructed without changing 
the basic Information found In the data that was not enhanced. This enhance- 
ment procedure along with a four-step check similar to that shown In figure 4 
was used to generate all of the remaining contour maps presented and to make 
certain that the enhancement process only clarified the contour plots and did 
0 not produce new features In the flow field. 


RESULTS AND DISCUSSIONS 

Introduction .- The data presentation Initially focuses on the use of con- 
tour maps of relative Mach number to show changes In flow conditions as the 
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rotor exit pressure was changed. Comparisons are made between the analysis 
code and the results obtained with the laser anemometer. Results are also 
presented In terms of overall pressure rise and turning as measured and pre- 
dicted. Plots of relative Mach number at various flow rates and percent spans 
are shown to give a more detailed picture of the flow In a quantitative sense. 
In each line plot the analysis Is compared to the LA data at comparable percent 
gap locations. Finally, some applications of the results which might have sig- 
nificance In the design of future compressor blading Is addressed. 

Mass flow rates measured In experimental facilities and calculated numeri- 
cally have a degree of uncertainty associated with them such that comparison 
of absolute magnitudes may not match well. To compensate for this fact, all 
of the comparisons are made at mass flow rates that have been nondlmenslon- 
allzed with the maximum flow rate measured from the experiment or calculated 
by the analysis code as appropriate. Although this approach was taken, the 
agreement between the predicted maximum flow and the experimental measurement 
was excellent for both the runs without and with boundary layer calculations. 
The maximum flow predicted without boundary layers was 1.2 percent higher than 
the measured flow and that predicted with boundary layers was only 0.06 percent 
higher. 


Relative Mach Number Contours 

Thirty percent span . - Figure 5 Is a composite showing the results ob- 
tained with the laser anemometer and analysis code at 30 percent span for the 
three flow conditions; max flow, peak efficiency, and the near-stall point. 

The shock locations are shown for each flow rate In order to provide a clearer 
picture of the flow In the passage. Locations of the wave systems and shocks 
were determined by Inspection of relative Mach number and relative flow angle 
data with respect to the streamwlse direction and designating the starting 
point of the flow deceleration as the shock front. This was done at as many 
pltchwlse locations as needed for the desired resolution. Experimental LA 
data must be Interpreted In this manner due to possible seed particle lag once 
a high velocity gradient region Is encountered. The numerical calculation 
actually produced shocks that were smeared over several grid points and thus 
the actual location of the shock Is between the beginning and ending of the 
deceleration. The grid spacing used herein was small, and the shocks were 
assumed to be located at the beginning of the deceleration. 

The LA data In figure 5 shows the progression from an oblique shock at 
max flow with' supersonic flow In the aft portion of the blading, to a near 
normal shock followed by a second shock at peak efficiency, and finally, a 
single normal shock ahead of the blade leading edge at the near stall point. 

The peak Mach number varies from 1.4 ahead of the shock at max flow and peak 
efficiency to 1.35 ahead of the shock at near stall. Also shown In the figure 
are the results obtained with the three-dimensional analysis code. The major 
difference between the experiment and analysis at max flow Is the location of 
the second shock at the exit of the passage. The peak Mach numbers upstream of 
the first shock are approximately equal. Although, the location of the second 
shock Is sensitive to backpressure at the max flow condition, a 7 percent In- 
crease In backpressure was required In order to move the computed shock to the 
approximate position of the measured second shock. This Is a large change In 
pressure, and It Is more probable that the difference In location Is a result 
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of the constant shape factor used In the boundary layer calculations. Also, 

LA data In the aft portion of the rotor for the max flow condition was only 
taken at the 40, 60, 90 and 100 percent chord locations so that good shock 
definition In this region was very difficult. The Mach number at the rotor 
exit Is less for the numerical results at this flow point than the measure- 
ments Indicate (an exit Mach number of about 0.95 for the data and 0.9 for the 
analysis) . 

The agreement at the peak efficiency point between the data Is much better 
than for the max flow case. The front shock Is accurately located by the code 
and there Is a rear shock located In approximately the same position as the 
measured rear shock. The second shock Is shown extending from the pressure 
surface to about midpitch. At points closer to the suction surface. Interpret- 
ing the location of the shock was more difficult since an acceleration of the 
flow or a constant Mach number region was not found In the streamwlse plots. 
However, Indications of a second shock were found near the suction surface by 
noting changes In the Mach number gradient In the streamwlse direction. This 
was true for both the analysis and the LA data. The computed Mach number In 
the exit region of the blade Is lower than the measured Mach number (0.8 com- 
pared to 0.9). Peak Mach number upstream of the shock Is 1.4 for both the 
data and the analysis. The agreement between analysis and data at the near 
stall point Is good In terms of shock location and peak Mach number. The com- 
puted Mach number at blade exit Is underpredicted (0.7 versus 0.8), as It was 
with the previous cases. 

Near tip and near hub sections . - Figure 6 shows the peak efficiency and 
near stall flow points at the 10 percent span location (fig. 6(a)) where the 
Inlet relative Mach number 1$ supersonic and at the 70 percent span location 
(fig. 6(b)) where the Inlet relative Mach number Is subsonic. Results are 
shown for both the experiment and the analysis. The peak efficiency LA data at 
10 percent span In figure 6(a) shows a front shock that Is nearly normal around 
the leading edge but Is oblique over most of the passage gap. Inspection of 
the streamwlse and pltchwlse LA data Indicated that a second shock might exist. 
However, the data was not conclusive and thus a second shock Is not Indicated 
In the figure. A second shock was predicted by the three-dimensional analysis 
code at the approximate location suggested by the LA data. The high values of 
relative Mach number after the normal portion of the shock have been Investi- 
gated by Straz1s?ar (ref. 5) arid appear to be due to the three-dimensional shock 
Inclination In the hub-to-shroud plane. 

At the near stall point the second shock has disappeared from the analysis 
and only a normal shock standing In front of the blade remains. A second shock 
system on the pressure surface around the leading edge Is Indicated by the LA 
data. This feature appears to be necessary to guide the high Incidence angle 
flow around the leading edge so that It can adjust to the blade surface angles. 
This assumption was supported by computing particle trajectories through this 
proposed system based on the wave angles and the measured Inlet Mach number. 
Final pressure surface flow angles after passing through the second shock sys- 
tem were within a few degrees of the surface values. This phenomenon Is not 
predicted by the analysis and may be due, at least In part, to tip region flow 
effects. More extensive analyses at other operating conditions and percent 
span locations with high Incidence angles Is required before a full explana- 
tion can be presented. Peak Mach numbers before the shock are, In both cases, 
approximately 0.1 higher for the analysis than for the data. Suction surface 
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Mach numbers are generally higher (0.1 to 0.2) In the analysis since the vis- 
cous effects are Incorporated only as a mass addition and not a momentum 
deficit. 

Figure 6(b) shows maps of constant relative Mach number for peak effi- 
ciency and near stall at the 70 percent span location for both the LA data and 
the three-dimensional analysis results. Flow field features between the two 
operating conditions are very similar and agreement between the LA data and the 
numerical results Is quite good although the analysis yields a slightly lower 
exit Mach number than experimentally verified. 


Overall Performance 

In figure 7 an attempt has been made to provide a clear picture of how 
the three-dimensional analysis code compares to the experimental data on a 
global basis that could be useful In the design phase of a new blade row. 

Figure 7(a) shows the static pressure rise across the rotor at the hub and tip 
as a function of . nondlmenslonal mass flow rate. As expected, the analysis 
without boundary layers predicted a much higher pressure rise at both hub and 
tip than was measured with conventional pressure Instrumentation. Excellent 
agreement between the measurements and the analysis with boundary layer effects 
Is experienced over the entire speed line from max flow to near stall. Un- 
stable operation of the rotor was experimentally observed at a mass flow ratio 
of about 0.89, while at a mass flow ratio of 0.906 the analysis was stable. 

The near-stall numerical solution was certainly aided by the restriction of a 
constant shape factor for the boundary layer calculations. Overall turning of 
the flow across the blade row at midpitch, given by rotor exit absolute tan- 
gential velocity divided by wheel speed Is shown versus percent span for the 
peak efficiency and the near stall point In figure 7(b). Agreement between LA 
data and analytical results Is within 30 percent at peak efficiency and less 
than 16 percent near-stall over the entire blade span. The pltchwlse averages 
of the flow turning (not shown) agree within 21 percent at peak efficiency and 
within 15 percent for near stall. The temperature ratios across the rotor 
predicted by the analysis and preliminary data from standard aerodynamic sur- 
vey probes (not shown) was compared at peak efficiency and found to agree with- 
in 2 percent over the entire span. This discrepancy between the LA results, 
analysis results, and standard. aerodynamic measurements Is expected to be re- 
solved by application of high frequency response pressure probes behind the 
rotor In future experiments. 


Relative Mach Number at MldPItch 

A more quantitative approach to comparison of the LA data and the analy- 
sis results Is to look at data along the rotor midpitch line. It Is precisely 
this type, of plot which has been used to determine the shock front locations 
for all of the relative Mach number contour maps presented thus far. The rela- 
tive Mach number along a midpitch line versus percent chord at 30 percent span 
for all three operating points Is shown In figure 8. The expansion of the flow 
after the first shock at max flow as predicted by the analysis code leads to a 
strong passage shock at approximately 70 percent chord. Some Indication of a 
second shock In the max flow LA data Is evident from the peak In Mach number 
at 10 percent chord, plateau from 30 to 40 percent chord and subsequent decel- 
eration from 40 to 60 percent chord. Agreement In the front shock location 
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between the analysis and experiment Is good but, as before, a low exit Mach 
number Is predicted by the analysis. The differences In the second shock loca- 
tion may be due In part to Insufficient LA survey locations In this region. 

Only one LA data point exists between the start of the computed reexpansion 
and the rotor exit. At peak efficiency the two-shock system that was observed 
In earlier contour maps Is again evident In the experimental data (shown by a 
plateau of Mach number from 20 to 30 percent chord) but not as evident In the 
analysis results. There Is however, a change In the gradient of computed Mach 
number at 40 percent chord Indicating the existence of the second shock. At 
percent gap locations closer to the pressure surface, the gradient change be- 
comes more pronounced until the existence of the second shock becomes evident. 
Shock location and peak Mach number before the shock agree well for both peak 
efficiency and near stall. 


Near-Blade Surface Comparisons 

Generally, an Important task for the blade row designer Is to tailor the 
blade shapes so that adverse gradients that produce large growths In the boun- 
dary layer are either reduced or controlled. It Is of great concern to deter- 
mine how well a particular analysis code predicts the surface velocities In 
order to determine how much faith the designer should place In the ability to 
tailor the blading for low losses. Blade loading In terms of differences In 
relative Mach numbers across the channel In the pltchwlse direction versus 
percent chord are shown In figure 9 for 10, 30 and 70 percent spans for both 
peak efficiency and near stall. The LA data Is taken at the first point off 
the blade surface where a sufficient number of measurements were available. 

This was generally close to 5 and 95 percent pitch and, therefore, the analysis 
results were selected at these same locations. For all three spanwlse loca- 
tions the agreement on loading levels Is reasonably good at peak efficiency. 

At near stall, however, the analysis yields much higher loadings at 10 and 
30 percent span but Is reasonable at 70 percent. The major difference In the 
loading level results from lower computed Mach numbers on the pressure surface 
than measured. Suction surface agreement Is generally good for both flow 
rates. For all percent span locations and flow points, the shock location on 
the suction surface Is well-predicted by the analysis. 

Another area which Is of some Interest to the turbomachinery designer Is 
the manner In which the flow Is turned through the blade row. Figure 10 shows 
the measured relative flow angle as a function of percent chord for both the 
suction and pressure surfaces. The results are shown for the 30 percent span 
Immersion location and for all three operating conditions. Also shown In the 
figure are the blade surface metal angles. The measured suction surface flow 
angles show that generally the flow follows the suction surface throughout the 
entire operating range. The analysis results of suction surface flow angle 
follow many of the same trends observed In the LA data. Flow deviations from 
the suction surface metal angle In the aft portions of the blade are observed 
both experimentally as well as analytically with the amount of deviation In- 
creasing from peak efficiency to near stall operation. The axial location at 
which the flow deviates from the blade changes from about 90 percent chord at 
peak efficiency to about 80 percent chord near-stall. The data at max flow Is 
too sparse to determine where the flow deviates from the blading. 

The pressure surface flow angle data In figure 10 appears not to be as 
well behaved as the suction surface data. Both analytical and experimental 
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data follow the blade surface over most of the chord, except In the leading 
edge region. Figure 10(c) Indicates that the pressure surface flow angle does 
not adjust and start to follow the pressure surface at near stall as quickly 
as It did at the other conditions. It Is not until about 30 percent chord 
that the relative flow angle matches the blade angle. The analysis also Indi- 
cated that the flow does not adjust to the blading Immediately at near-stall* 
however It adjusted more rapidly than the LA data Indicates. 


Analysis Results In the Meridional Plane 

Figure 11 shows the nature of the shock system and how the numerical 
result changes from peak efficiency to near stall. In this suction surface 
hubto-shroud view the shock Inclination angle changes dramatically at radial 
positions between the tip and 30 percent span while below 30 percent span, the 
Inclination remains constant. The shock location changes from the aft portion 
of the blading near the tip to the front part of the blading near the hub. 
Overall, the Mach numbers in front of the shock Increase as the operating point 
changes from peak efficiency to near stall. Since It Is difficult to predict 
the point where stall will occur In the analysis, due to the use of a constant 
shape factor, one approach Is to consider where the shock structure Is strong 
enough to separate the boundary layer. The static pressure rise across a nor- 
mal shock at a Mach number of 1.3 is generally considered sufficient to produce 
boundary layer separation. At the near stall flow, the Mach number contour at 
1.3 penetrates to about 90 percent span from the tip and the Inclination of 
the shock In the meridional plane decreases from that at peak efficiency. The 
shock ,1s almost normal to the flow at near stall, and thus, the pressure rise 
across It Is approaching that of a normal shock. The shock location has moved 
forward and occurs In front of the rotor throat. This combination could pro- 
duce global separation on the blade and subsequent stall of the rotor. Infor- 
mation of this kind could be used for a new design to warn the designer of a 
potential stall condition. 


CONCLUDING REMARKS 

The comparisons between analytical and experimental data for a transonic, 
low-aspect ratio, axial-flow fan rotor have, In general, shown good agreement 
for all of the operating conditions presented. The overall trends In the rela- 
tive Mach number contour maps Indicate the very good shock capturing ability 
of the analysis code at all rotor span locations. The results presented have 
demonstrated the extent to which an Invlscld Euler solver with a very simple 
boundary layer correction scheme was able to predict the complex three- 
dimensional flow field of this transonic rotor flow field. The code predicted 
the maximum flow for the rotor extremely well. Results were also shown that 
Indicate that the code may be useful In analysis of a new blade row design to 
Insure that desired stall margin Is attained. The chordwlse plots of relative 
Mach number and relative flow angle reveal some areas In which analysis and LA 
data deviate. Some of these conditions may be traced to the boundary layer 
model used. 

A three-dimensional analysis technique, such as the one used herein, would 
probably not be used In the preliminary design of turbomachinery blading, how- 
ever It may well be used to check a final configuration for design and off- 
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design operation. This capability to compute the off-design flow field Is 
extremely Important and has shown to yield reasonably good quantitative results 
for the present rotor configuration. The analysis code did not predict the 
multiple shock system on the pressure surface required to turn the flow paral- 
lel to the blade. This may have been due to the numerical smoothing used In 
the code and should be Investigated further. Results concerning the Intra- 
blade flow turning and work distribution should be useful Information for com- 
pressor design systems. 

As mentioned, the goal of the current NASA/ARMY research program Is cen- 
tered around Increasing the fundamental understanding of the three-dimensional 
flow field of compressor blading. The material presented has been one such 
attempt to Investigate the flow field of a particular transonic fan rotor, and 
as such, the results shown apply only for this case. Further work dealing 
with different types of blading Is needed to determine how much of what was 
observed experimentally and computed numerically for this Isolated rotor. Is 
generic to this class of turbomachine and how much was rotor specific. As 
more machines are Investigated , questions such as blade row Interactions and 
multistage matching should also be addressed. 
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(a) Contour plot using "raw" data. (b) Conour plot with all "zero" data removed. 



tc) Contour plot with data interpolated to increase td) Contour plot with final level of smoothing, 
the number of axial locations. 

Figure 4. - Contour plots of relative Mach number at 30% span from the tip at the peak efficiency 
flow rate showing the effect of the enhancement process on the LA survey data. (Note: blades 
not shown. I 
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Figure 5. - Contour plots of relative Mach number at 30% span from the tip for max flow, peak 
efficiency, and near stall flow points: comparison of LA survey results with analysis results. 
— Shock location. 
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(a) 10% span. 

Figure 6. - Contour plots of relative Mach number at 10% span from the tip for the peak efficiency 
and near stall flow rates: comparison of LA survey data with analysis results — shock location. 
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(b) 70% span. 

Figure 6. - Concluded. 
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Figure 8. - Variation in relative Mach number with percent chord at mid-pitch 
for three flow points at 30*percent span from the tip. 
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(a) 10-Percent span from tip, (b) 10-Percent span from tip, near 

peak efficiency. stall. 



(c) 30-Percent span from tip, (d) 30-Percent span from tip, near 

peak efficiency. stall. 
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(e) 70-Percent span from tip, (f) 70-Percent span from tip, near 

peak efficiency. stall. 

Figure 9. - Variation in relative Mach number with percent chord at three 
spanwise locations for peak efficiency and near stall flows. 
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Figure 10. - Variation in relative flow angle with percent chord along the speed line for experiment and analysis 
at 30-percent span from the tip. 




Figure 11. - Contours of suction surface relative Mach number projected on the meridional hub-to-shroud plane, 
analysis results only. 
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